Azimuthally sensitive correlations in nucleus-nucleus collisions 
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We define a set of model-independent observables which generally characterize the azimuthal 
dependence of two-particle correlations in nucleus-nucleus collisions. We explain how they can 
be analyzed, and show to what extent such analyses are model dependent. We discuss specific 
applications to the anisotropic flow of decaying particles, azimuthally sensitive HBT, and correlations 
between particles at large transverse momentum. A quantitative prediction is made for jet quenching 
with respect to the reaction plane. 
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PACS numbers: 25.75.Gz, 25.75.Ld, 25.75.Nq, 12.38.Mh 



I. INTRODUCTION 

In non-central nucleus- nucleus collisions, azimuthal an- 
gles of outgoing particles are generally correlated with 
the direction of the impact parameter. This phe- 
nomenon, called "anisotropic flow," has been known 
for 20 years Q, and has raised particular interest at 
the Brookhaven Relativistic Heavy Ion Collider (RHIC) 
where it is thought to provide unique evidence for quark- 
gluon plasma (QGP) formation 0|. 

Most often, one studies the azimuthal dependence of 
single-particle production [3(|. Here, we would like to dis- 
cuss the azimuthal dependence of two-particle correla- 
tions. This is of interest in various situations: 

• Anisotropic flow of short-lived particles: the flow 
of unstable particles (for instance A baryons) is 
studied through their decay products. One must 
first identify a correlation between daughter par- 
ticles, typically through an invariant mass plot; 
then study how this correlation depends on the az- 
imuthal angle of the decaying particle 0, IE IE IE IE 



9j. 



• Azimuthally sensitive two-particle interferometry: 
Bose-Einstein correlations between identical parti- 
cles are commonly used to measure the size and 
shape of the emitting source 01 ■ I n non-central 
collisions, the source projection on the transverse 
plane is no longer circular 

El 

, and this can directly 
be seen in Hanbury-Brown Twiss (HBT) studies of 
two-particle correlations, as already observed at the 
Brookhaven Alternating Gradient Synchrotron [l^ 
and at RHIC HE]. 

• Jet quenching with respect to the reaction plane: 
the energy loss of hard partons traversing a decon- 
fined medium 113 is a crucial signature of QGP 
formation at RHIC |16|. In particular, it results 



in a modification of the pattern of azimuthal cor- 
relations between high-py hadrons, compared to 
pp collisions: the back-to-back correlation is sup- 
pressed ElEj. In a non-central collision, the aver- 
age length of matter traversed by a parton depends 
on its azimuth |l9ll20| . which results in azimuthally 
dependent two-particle correlations [2lL 12^ . 
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In this paper, we give for the first time a unified pre- 
sentation of these phenomena, which have so far been 
discussed separately. In all cases, analyzing azimuthally 
dependent correlations involves two distinct operations: 
1) Measuring the distribution of a pair of particles with 
respect to the reaction plane; 2) Isolating the "true" cor- 
relation from the uncorrelated part. Both issues can be 
discussed independently, on a fairly general footing. 

The first operation is discussed in Sec. |H] where the 
observables associated with two-particle anisotropic flow 
are defined. These observables are model-independent 
and can in principle be measured accurately. In particu- 
lar, it will clearly appear that any method used to mea- 
sure the single-particle anisotropic flow can also be used 
to analyze azimuthally sensitive correlations, modulo mi- 
nor modifications. In that view, we recall in Appendix 1X1 
the main features of existing methods for analyzing one- 
particle flow, and we introduce the changes necessary to 
measure pair flow. While existing methods all require 
to estimate the reaction plane on an event-by-event ba- 
sis HE Hi (see for instance Refs. E|l for A flow, HJ 
for azimuthally sensitive HBT, and [27| for correlations 
between high- momentum particles), this step is by no 
means necessary with the procedure we suggest. This 
opens the possibility to apply the improved methods of 
flow analysis recently devised in Refs. [28L 12^. I30L l3lll33| 
and to resolve an inconsistency of present analyses: on 
the one hand, one studies a correlation (between decay 
products, due to quantum statistics, from jet fragmenta- 
tion) which is essentially a "nonflow" correlation; on the 
other hand, one uses the event-plane method which relies 
on the assumption that all correlations between particles 
are due to flow |23| . 

The second operation is discussed in Sec. IIIII Un- 
like the first one, it will be shown to be always model- 
dependent. Several specific applications are discussed in 
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Sec. lIVI together with predictions regarding the pair-flow 
coefficients. Our results are summarized in Sec. ]V\ 



II. OBSERVABLES FOR TWO-PARTICLE 
ANISOTROPIC FLOW 

We first recall definitions for single-particle distribu- 
tions. For particles of a given type in a given rapidity 
(y) and transverse momentum (j>t) window, the prob- 
ability distribution of the azimuthal angle 4> (measured 
with respect to a fixed direction in the laboratory) reads 
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FIG. 1: Illustration of the various azimuthal angles <f>i, (j>2, 
pa ir, A(j>, with x = §. 



where is the (unknown) azimuth of the reaction plane 
(impact parameter) in the laboratory frame. 

The Fourier coefficients 33] in this expansion are given 
by v n = (e - " 1 ^ - **)), where angular brackets denote an 
average over particles and events. Given the normaliza- 
tion choice in Eq. (JJJ, vq — 1. Since p(<f> — $jj) is real, 
V- n = (v n )*, where * denotes the complex conjugate. If, 
in addition, the system is symmetric with respect to the 
reaction plane [—(</> — &r) is equivalent to <j> — Qr], as 
in a collision between spherical (although not necessarily 
identical) nuclei when parity is conserved, Eq. (JIJ reads 



1 
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with v n = (cosn((f) — $>r)), i.e., v n is real. 

The Fourier coefficients v n have by now become famil- 
iar in the study of anisotropic flow in (ultra)relativistic 
heavy-ion collisions. Nevertheless, it is instructive to re- 
call why they are the proper tools to parameterize az- 
imuthal anisotropics. The key feature is that even though 
the reaction plane <&r is unknown on an event-by-event 
basis, the first Fourier coefficients v n can be accurately 
reconstructed from a statistical analysis of azimuthal cor- 
relations between outgoing particles (see Appendix IA II 
for a review of the methods for analyzing single-particle 
flow). However, the higher the value of n, the larger 
the uncertainty on v n 34]. Therefore the probability 
p((j) — $/f) at a specific azimuth cannot be measured in 
practice. Furthermore, since v n is defined as an aver- 
age, it is also easier to compute in theoretical studies — in 
particular, in Monte Carlo models — than the probability 
distribution itself. 

The above definitions can readily be generalized to the 
distribution of particle pairs with respect to the reaction 
plane. A pair of particles of given species is character- 
ized by 6 kinematic variables pn, 2/i, 4>i,PT2> 2/2, 4>2- It 
is convenient to combine <pi and </> 2 into the relative an- 
gle A(j) = 4>2 — 4>\ (or any similar observable that does 
not depend on the overall orientation of the pair in the 



transverse plane, as e.g. the invariant mass) and a "pair 
angle" 



X(j>\ + (1 - x)4> 2 , 



(3) 



where < x < 1. One can restrict </> pa ir and A<p to the 
ranges — ir < pa i r < ir and — ir < A<p < ir. If x = \, 
^pair is the mean angle. The choice of x depends on the 
problem under study: most often, one chooses for </> pa ir 
the azimuthal angle of the total transverse momentum 
Pti +PT2 ( see Sec. II V Al and Sec. HVB|) ; in studies of az- 
imuthal correlations between high-momentum particles, 
x = 1 is a more common choice (see Sec. IIV Cjl . 

Consider now a sample of pairs of particles in some 
range of Pt%, PT2, 2/i, 2/2, Acf>. To study the probability 
distribution of the pair angle pa i r within this sample, we 
write its probability distribution in a way analogous to 
Eq. 0}: 



p(0 P air - $r) = 



+ 00 



-*h) 



(4) 



As the usual v n 's, the "pair- flow" coefficients v£ an are 
defined by «P air = ( e - in (^-*«)), with the normaliza- 
tion fg air = 1. Since the probability distribution is real- 
valued, the coefficients also satisfy the property 



pair 



« ai T 



(5) 



But unlike the single-particle flow v n , the pair-flow coeffi- 
cient v P air is in general not a real number. The underlying 
reason, exemplified in Fig. ^ is that the transformation 
0pair — &R — > — (0 P air — $r) for a constant A(j> is not a 
symmetry of the system. 1 As a consequence, sine terms 



The actual symmetry is under the simultaneous transformation 
</> P air - $i? — » -(^pair - — > —A<j>. Its consequences for 
the coefficients «£ alr are discussed in Appendix ITU 
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are also present in the real form of the Fourier expansion, 
and Eq. (J2J is replaced by 

1 ( + °° 

K^pair - *r) = Y 1 + 2 5] W cos n (<f>P™ - $ «) 
\ n=l 

+ <f smn(^ pair -$ fi )]\ (6) 

where the real coefficients v P ^ n = (cos n(0 pa i r — $^)) and 
v s?n = (sinn(0 P air — are related to the complex v^ alT 
by the relation v^ a " = w P ^ r - fu P ^ r . 

The existence of such sine terms was already noted 
in Ref. |2|j in the context of azimuthally sensitive HBT 
studies. Note that it does not imply parity violation, as it 
would in the case of single-particle flow [35| . The physical 
meaning of these additional terms will be illustrated in 
Sec. IIVI in various physical situations. In particular, we 
shall show that they may yield insight on the mechanism 
responsible for the deficit in high-p^ particles. 

It may be interesting to note that the pair-flow coeffi- 
cients vP air (A(f>) , when viewed as functions of the relative 
angle A</>, have peculiar properties which are derived in 
Appendix [BJ Checking that measured values of the co- 
efficients possess these properties then provides a way to 
evaluate the errors affecting the measurement. 

In experimental analyses, any method that can be used 
to measure the single-particle flow v n can be applied to 
extract the cosine terms v^ a " , without any modification: 
one simply considers the pair as a single particle with 
azimuthal angle pa ir- The generalizations required in 
order to extract the sine terms v ps ^ r are quite straightfor- 
ward. They are summarized in Appendix IA 21 for various 
methods of flow analysis. In particular, some of these 
methods can safely correct for nonflow effects, others for 
acceptance anisotropies. This extends the possibility of 
studying azimuthally-dependent correlations to detectors 
with partial azimuthal coverage. 

To conclude this Section, let us emphasize that 
the new characterization of azimuthally-sensitive two- 
particle correlations which we propose, with pair-flow 
Fourier coefficients, represents in our view an improve- 
ment over previous parameterizations, in the same way 
as v n is an improvement over older observables for 
anisotropic flow. The reason is simply that v? a ™ an d w p „ r 
are model-independent and detector-independent observ- 
ables. 



III. ISOLATING THE CORRELATED PART 

Subtracting the "trivial" uncorrelated part in order to 
isolate the "true" correlation is far from trivial. In this 
Section, we discuss this issue in as simple and general a 
way as possible. For the sake of simplicity, we start with 
the case when there is no anisotropic flow. In Sec. IIII Al 
we explain why the subtraction always involves some de- 
gree of arbitrariness, most often in the form of an ar- 



bitrary constant. Although this is to some degree well- 
known, at least to those who actually perform correlation 
analyses, we think it is worth recalling, since the litera- 
ture on the subject is rather confusing. In Sec. IIII Bl 
we recall the various ways of normalizing the correlation, 
depending on the observable under study. 

In practice, however, anisotropic flow is most often 
present, and makes the background subtraction more dif- 
ficult in heavy ion collisions than in elementary collisions. 
Note that this applies to all correlation analyses, not only 
to the azimuthally-dependent ones: the correlation of sin- 
gle particles with the reaction plane induces a correlation 
between them, which must be always subtracted, at least 
in principle, in order to isolate other effects (see for in- 
stance |3(j). In Sec. IIII CI we explain how this can be 
done, and show that this subtraction implies further ap- 
proximations. 



A. A model-dependent issue 

In a given event, let N± and denote the numbers 
of particles in two phase-space bins, (pn , J/i j 0i ) an d 
(PT2, 2/2, 4>2)- To simplify the discussion, we assume that 
the two bins are separated. If they overlap, one need 
only replace N\Na by the number of pairs in what fol- 
lows. The simplest definition of the correlation between 
the two bins is 

C = (NtN 2 ) - (JVi) (N a ) , (7) 

where angular brackets mean an average over many 
events. 2 

Such a definition is not satisfactory in practice because 
the sample of events used in the analysis always contains 
events with different centralities: in particular, the to- 
tal multiplicity may have sizeable fluctuations within the 
sample of events considered, and these fluctuations alone 
induce a correlation between any two phase-space bins. 
This correlation is of a rather trivial nature, but it may 
well overwhelm the interesting ones |37|. A simple way 
out of this problem would be to normalize the two terms 
in Eq. J7J) by the total number of pairs of correlated and 
uncorrelated particles, respectively, that is, to define in- 
stead the correlation as 

C = {NiN 2 ) - {N{ " n ~ 1)} (iVr) (N 2 ) , (8) 

where N is the number of particles in a large enough 
phase-space bin (typically, the total number of charged 
pions seen in the detector for interferometry analyses). 
This definition is also unsatisfactory for two reasons: 



2 In other terms, (N\), {N2) and {N1N2) are the one- and two- 
particle inclusive cross-sections, divided by the total inelastic 
nucleus-nucleus cross-section. 
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first, it obviously introduces some degree of arbitrari- 
ness in the definition of the correlation, depending on 
the choice of the phase space for N; second, part of the 
fluctuations in N may be meaningful for the correlation 
analysis, as in the case of Bose-Einstein correlations [38l| . 
so that there is no point in subtracting them. 

In actual analyses, the correlation is rather defined as 

C=(N 1 N 2 )-c(N 1 )(N 2 ), (9) 

where c is some free coefficient. This coefficient is kept 
constant throughout the correlation analysis (which typ- 
ically involves varying the invariant mass, the relative 
momentum, or the relative azimuth between the two par- 
ticles) . It is then fitted in such a way that the correlation 
C vanishes when it is expected to: at large relative mo- 
mentum in HBT analyses |38| . in some range of A(j> in 
correlations between high-p^ particles |171 Il8j . In addi- 
tion, as mentioned above, the background subtraction is 
often complicated by the existence of anisotropic flow, as 
we shall discuss in Sec. 1111 CI 

B. Normalizations 

There are essentially three ways of normalizing yields 
of particle pairs in azimuthally-independent analyses, de- 
pending on the observable under study. 

1. One simply computes the average number of pairs 
per event, {N\N 2 ). For instance, in order to mea- 
sure A production, one plots the number of (p, 7r _ ) 
pairs per event as a function of the invariant mass 
M of the pair. The number of pairs in the peak 
around the A mass gives the yield of A baryons, 
modulo acceptance corrections. 

2. One divides the number of pairs (N1N2) by the 
number of uncorrelated pairs, c{Ni)(N 2 ). This is 
the standard observable for Bose-Einstein correla- 
tions, where the ratio varies ideally between 2 and 
1 as the relative momentum of the pair increases. 

3. The third, intermediate choice is to divide the num- 
ber of pairs (Ni N 2 ) by the number of "trigger par- 
ticles" (Ni) [T3I . After subtraction of the uncorre- 
lated part c(N 2 ), one thus obtains the mean num- 
ber of particles N 2 correlated with a trigger parti- 
cle, which is independent of the system size (i.e., the 
same for a nucleus-nucleus and for a proton-proton 
collision) if there is no final-state interaction. 

C. Subtracting the correlation due to flow 

When the particles in the pair are individually corre- 
lated with the reaction plane this induces a trivial 
correlation between them, which must also be subtracted. 

This subtraction is easy in principle: one simply re- 
peats the operations of Sees. IIII Al and IIII Bl for a fixed 



orientation of the reaction plane Then, the following 
substitutions hold 

(JViJVa) {N 1 N 2 )(2i:) P ((b pail -^ R ) 
(JVi) (JVi)(27r)pi(^i-$ fl ) 
(AT 2 ) _> (N 2 )(2Tr)p 2 {<f> 2 -<Z> R ). (10) 

In these equations, (NiN 2 ), (Ni) and {N 2 ) denote quan- 
tities averaged over $r; p(0 pa ir — $i?) is the distribution 
of the pair angle, defined in Eq. iJBJ, and pi(<j>i — §r), 
Viifyi — <&r) denote the single-particle azimuthal distri- 
butions of each particle, defined as in Eq. 

Once the azimuthal distributions of pairs and single 
particles with respect to the true reaction plane have 
been properly reconstructed, extracting the correlation, 
and its azimuthal dependence, is straightforward. 

Strictly speaking, however, the one- and two-particle 
probabilities p(<f) — Qr) and p(0 pa ir — &r) at a specific 
azimuth relative to the reaction plane cannot be recon- 
structed. As already mentioned in Sec. [HJ only the first 
few Fourier coefficients v n or «P air can be reconstructed, 
due to larger absolute uncertainties on higher-order co- 
efficients. On the other hand, the Fourier coefficients of 
a smooth function of 4> — ®R are expected to decrease 
quickly as the order increases (this expectation is sup- 
ported by recent experimental single-particle flow data 
[33), so that one can reasonably truncate the series, 
keeping only the measured coefficients. This truncation is 
always required in order to estimate the correlation from 
anisotropic flow [rR l27j. At RHIC, for instance, the er- 
ror on the azimuthal distribution at midrapidity is likely 
to be dominated by the error on the fourth harmonic 1)4, 
and one can take ^5v4 as the error on p((f> — &r). 

IV. APPLICATIONS 

We shall now discuss specific applications, with em- 
phasis on the details of the experimental procedure. 

A. Anisotropic flow of short-lived particles 

Let us begin with the measurement of the anisotropic 
flow of particles that are seen through their decay prod- 
ucts, such as A ^ptt" H0EI3, 71-0 77; V -> 77 0: 
Kg — > 7T + 7r - |(|, [ij. We shall illustrate the recipe by 
discussing the flow of A baryons. 

For each event, one sorts (n~,p) pairs into bins of 
invariant masses M. The first step is then to analyze 
the total A yield: following the standard procedure, one 
counts the number of pairs in each invariant-mass bin, ir- 
respective of the pair azimuth; let -/V pa i rs (M) denote this 
number. One then separates this distribution into an un- 
correlated part [the background Nt,(M)] and a correlated 
part [the peak Na(M), centered around the expected A 
mass] : 

N pails {M) = N b (M) + N A (M). (11) 



5 



The integral of the correlated part N\(M) over M is the 
A yield. 

In most cases, the peak is well above the background: 
to perform the above decomposition, one need not go 
through the whole procedure of the previous section: in- 
stead, one simply assumes that the background Nb(M) 
is a smooth function of M 25]. Please note that the 
anisotropic flow of n~ and p correlates their azimuthal 
angles, and therefore distorts the background. However, 
the distorted background remains smooth, so that this 
effect need not be taken into account. 

Next, one defines the azimuthal angle of the pair, pa i r , 
as the azimuthal angle of the total transverse momentum 
PTi + PT21 and one analyzes the pair flow coefficients 
v^"(M) and v^ r (M) in each bin. One then performs 
a decomposition similar to Eq. for the azimuthally 
dependent part of the pair yield: 

iV pairs (M) v c , n (M) = N b (M) »W (M) + N A (M) « C A „ 
iV pail ,(M) v s>n {M) = N b (M) vi b l(M) + N A (M) v^ n . 

'(12) 

This decomposition is performed assuming that the back- 
ground components N b (M)v£l(M), N b (M)v { s b) n (M) are 
smooth functions of M. In this particular case, sym- 
metry with respect to the reaction plane for A particles 
implies v^ n — 0, except for experimental biases and fluc- 
tuations. This identity can be used in order to check the 
accuracy of the experimental procedure, as in the case 
of single-particle flow j24[. If the background consists of 

uncorrelated particles, one also has ««?n(M) = 0. 

In most analyses so far, the decomposition between the 
background and the peak is performed independently for 
several bins (typically, 20) in pair — [S El EH) where 
tyji is an estimate of the reaction plane. With the above 
procedure, the decomposition is only performed twice, in 
Eqs. (jTTj) and flU) . 

When the peak/background ratio is low, finally, mixed 
events can be used to define the background [5j. However, 
the above-mentioned distortion of the background due to 
anisotropic flow must then be taken into account, as we 
shall see in more detail in Sec. IIVBI 



B. Quantum correlations 

Azimuthally dependent Bose-Einstein (or, more gener- 
ally, short-range) correlations are analyzed in two steps. 
The first step is to perform a Fourier expansion of the pair 
yield with respect to the reaction plane, for each relative 
momentum q |HJ. As explained in Sec.[H] any method 
of flow analysis can be used to extract the Fourier coef- 
ficients w c ,n(q) and w Si „(q). Even with the event-plane 
method, no binning injA pa jr — is required, in contrast 
to present analyses [Tj, Il3| . Once the coefficients are 
known, one easily builds the distribution of pairs relative 
to the reaction plane (up to truncation issue mentioned 
in Sec lHTCTl . 

Next comes the difficult part of the analysis: one must 
divide the number of pairs per event by the number of 
uncorrelated pairs, as explained in Sec. lIIIBl For a fixed 
orientation of the reaction plane <I>/j, this number de- 
pends on $ fl (see Sec. lIIICTl : 

AUcor($i?.) = 27r(JVi)p 1 (^i-$ fl ) 

X27T(iV 2 }p 2 (</» 2 -$ fl ). (13) 

The ^-independent part, (Ni) (A" 2 ), can be obtained us- 
ing a standard mixed-event analysis. The <I>fl-dependent 
part, however, involves the (first) flow coefficients v n of 
both particles in the pair. 

To avoid this complication, the procedure suggested 
in Ref. [2(| is to used mixed events with aligned event 
planes. This procedure, however, is only approximate, 
because one mixes events with different reaction planes, 
although the estimated planes are the same. To be 
specific, let us compare in a simple case the distribu- 
tion of uncorrelated pairs following the exact procedure, 
Eq. I|13(l . and using mixed events with aligned event 
planes. To simplify the calculation, we assume that only 
elliptic flow v 2 is present, and that it has the same value 
for both particles in the pair; we further assume that the 
pair angle is the mean angle, x — ^ in Eq. ©■ Then the 
exact result is 



N ZTft^ = 1 + 2v 2 cos2A0 + Av 2 cos2A0 cos2(0 pair - + 2v\ cos4(0 pair - (14) 

{1\2) 

This is to be compared with the result obtained following the method of Ref. [2(j : 

^Tw^ = l+2v l (cos2Avl/ J ,) 2 cos2A0 + 4^ 2 cos2A0 cos2(0 pair - ^ R ) + 2v 2 2 ^^f^ \ cos4(^ pair -$ fl ), (15) 
{I\2) (C0S4AWfl) 



where A^j? = — is the difference between the estimated event plane and the true reaction plane. As 



expected, both results coincide when A^r = 0. Quite 
remarkably, the mixed-event method is correct to lead- 
ing order in vi even when A^^j 7^ 0. However, it misses 
the coefficients of order v\. Typical values of the cor- 
rection factors for the STAR experiment at RHIC are 
(cos2A<Fr) 2 ~ 0.6 and (cos2A1^) 2 / (cos4A* fl ) ~ 1.3. 

Besides the systematic uncertainty we just discussed, 
the price to pay for aligned mixed events is that one 
must essentially perform the whole correlation analysis 
for fixed values of both the pair angle pa ir and the es- 
timated reaction plane ^r. We suggest instead the fol- 
lowing method: 

1. place particle pairs in bins according to their ra- 
pidities yx, y2, total transverse momentum K = 
Pti + PT2: an( i relative momentum q; 

2. in each such bin, build the correlation function 
C(q) as in the standard, azimuthally-insensitive 
HBT analysis; 

3. reconstruct the azimuthal distributions of pairs, 
p(^ P air - and of single particles, p^fa - $ R ) 
and P2{4>2 — *&r)) with respect to the actual reac- 
tion plane; that is, measure the first Fourier coeffi- 
cients wP air and v n for each particle in the pair (at 
RHIC energies, measuring the second and fourth 
harmonic, «2 and U4, should be enough to guaran- 
tee that the distributions are reasonably well recon- 
structed); 

4. with the help of the substitution Eq. 1)100. build the 
azimuthal dependence of the correlation function. 

One then eventually extracts azimuthally-dependent 
HBT radii using standard techniques (in particular, in- 
cluding correction for Coulomb effects) which are beyond 
the scope of this paper. 



C. Two-particle azimuthal correlations 

Two-particle azimuthal correlations at large transverse 
momentum are under intense investigation in ultrarela- 
tivistic nucleus-nucleus collisions, since it has been re- 
alized that they yield direct evidence for hard scatter- 
ing j3||. In that case, one correlates a high-pr particle, 
the "trigger" particle, hereafter labeled 1, with a lower- 
Pt particle, hereafter labeled 2. We assume for simplic- 
ity that particles 1 and 2 belong to separate px intervals. 
This is not the case in the STAR analysis [l7|. where 
particle 2 can be any particle with momentum lower than 
Pt\ above some cut. This difference, however, is not cru- 
cial for the following discussion. 

The following quantities must be measured: the aver- 
age number of pairs per event, as a function of the relative 
angle A.<f>, (iVp a j rs (A^>)) (in practice, pairs are naturally 
sorted into equal-size bins of A(f>) , and the average num- 
bers of particles per event (Ni) and (N?). 
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FIG. 2: Illustration of the prediction made in Eq. 1171 . The 
long arrows represent the momenta of trigger particles, while 
the shorter arrows represent the momenta of associated par- 
ticles. If the modification of the correlation is due to the 
quenching of the associated particle, it must be unchanged 
under the transformation <f>2 — &r — > — (<f>2 — &r)- 

In studying the azimuthal dependence of the correla- 
tion, a natural choice is to take the azimuthal angle of the 
trigger particle <pi as the pair angle pa i r , he., one chooses 
x = 1 in Eq. (J3J. One needs to reconstruct the azimuthal 
distribution of pairs (for a given A</> bin), p A ^{(f)i — &r), 
and the azimuthal distributions of both trigger and as- 
sociated particles, Pi(4>i — &r) and ^2(^2 — $r)- One 
may then reconstruct the whole correlation function for 
a fixed value of <f>\ — $>r. In particular, the correlation 
functions for the specific values <pi = (in plane) and 
4>i = $r + tt/2 (out of plane) are given by 

c ,„ (A ^ )= W^))Ao) _ 2;r 

(iVi) _pi(0) 

where c is a constant close to unity, as explained in 
Sec. IIII Al It is independent of A0 and 4>i — 

Let us briefly compare the above procedure with the 
one suggested by Bielcikova et al. [23: these authors 
show how to analyze correlations in- and out of an event 
plane, which is not exactly the reaction plane. Since the 
event-plane resolution is a detector-dependent quantity, 
this prevents quantitative comparisons between different 
experiments. In addition, the algebra to subtract the 
uncorrelated part is much simpler with our method. 

The standard interpretation of the modification of the 
correlation function in nucleus-nucleus collisions, com- 
pared to proton-proton collisions 01 > * s that the asso- 
ciated parton loses energy on its way through nuclear 
matter. If this interpretation is correct, then for a given 
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A<f> the number of pairs per trigger particle depends only 
on the path followed by the associated particle. Symme- 
try with respect to the reaction plane implies that it is 
unchanged if fa — &r is changed into its opposite. As 
illustrated in Fig. [2J this symmetry is by no means triv- 
ial since the path followed by the trigger particle is now 
different. This gives us, for arbitrary fa — &r and Afa 
the prediction 

p A4> (fa - $7? - Afa = P A H-4>2+^R~ Afa {m 
Pi (02 -*K- A0) Pl(-fa + ®R- Afa ' 

If the only nonvanishing Fourier harmonic in the single- 
particle and pair azimuthal distributions is i>2, a simple 
calculation shows that the previous identity is equivalent 
to 

<f (A0) = (<f(A0) - tan2A0, (18) 

where i»2 is the elliptic flow for the trigger particle, and 
the pair-flow coefficients w^ 2 J and a 2 ir have been defined 
in Sec. [HI This prediction is consistent with the general 
symmetry property (|B3|) . 

V. SUMMARY AND PERSPECTIVES 

We have introduced novel, model-independent observ- 
ables that describe the dependence in azimuth of two- 
particle correlations in heavy-ion collisions. These ob- 
servables, namely the coefficients v^ T and v^™ 1 in the 
Fourier expansion of the azimuthal distribution JSJ of the 
pair-angle </> pa ir that characterizes (together with the rel- 
ative azimuth) particle pairs, generalize in a natural way 
the Fourier coefficients v n for single-particle anisotropic 
flow. As the latter, the pair-flow coefficients can easily be 
measured in experiments, using any "usual" method of 
flow analysis (modulo minor modifications for the mea- 
surement of the sine terms, v^ a "): event-plane, two- 
particle correlations, cumulants, Lee- Yang zeroes can 
equally be applied. We however recommend the last two, 
when possible, in order to disentangle flow from non-flow 
effects. 

A main point of this paper is that these observables 
should replace, in future analyses, quantities that are 
defined for a given azimuth relative to the event plane. 
Much in the same way, the Fourier coefficients v n have 
now replaced earlier observables such as the "flow an- 
gle" , the "squeeze-out ratio" in most, if not all analy- 
ses of single-particle flow. The reason is that the event 
plane is not exactly the reaction plane, and the disper- 
sion varies from one experiment to the other. Therefore, 
observables defined with respect to an event plane only 
yield a qualitative information. The v£ air , on the other 
hand, allow studies of azimuthally-sensitive correlations 
to enter the quantitative era. 

In a second part, we have briefly shown how to re- 
late our observables to physical quantities of interest in 



three different cases: anisotropic flow of decaying parti- 
cles, interferometry, and azimuthal correlations of high- 
momentum particles. It is important to stress that, un- 
like the measurement of the pair-flow coefficients v?™ 1 
and uj^f, this second step does depend on the under- 
lying physical picture. This model-dependence leads to 
some arbitrariness, which in practice takes the form of the 
introduction of a normalization constant and (for HBT 
and high-py-particle studies) a necessary truncation of 
the Fourier expansion of the single-particle distribution. 

A striking difference with usual studies of anisotropic 
flow is the general occurrence of sine terms in the Fourier 
series expansion. The relevance of such sine terms was 
already discussed in the context of single-particle flow 
[24l |33 and azimuthally sensitive HBT (2|j. Here, we 
have shown for the first time that in the case of correla- 
tions between high-p^ particles, jet quenching would re- 
sult in a specific value for the sine term, given by Eq. 11811 . 

Azimuthally-sensitive correlations are among the sub- 
tlest analyses in our field, and have already given valu- 
able insight on the physics of high-energy nuclear colli- 
sions. We hope that the observables and methods in- 
troduced in this paper will help to improve future analy- 
ses. Thanks to the high statistics now available at RHIC, 
or that one can anticipate at the CERN Large Hadron 
Collider (LHC), new measurements will become possible, 
for instance the azimuthal dependence of non-identical- 
particle interferometry or the anisotropic flow of various 
"new" particle types, while probing new regions of phase 
space. 
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APPENDIX A: METHODS FOR ANALYZING 
SINGLE- AND TWO-PARTICLE FLOW 

As stated in Sec. |Hj the measurement of the pair-flow 
coefficients v? a " an d v^ a " involves the same methods of 
analysis as for single-particle flow coefficients v n (mod- 
ulo a small modification when measuring v^ a "). This 
prompts us to recall briefly the various methods that have 
been proposed in the literature, indicating the modifica- 
tion necessary to measure the sine coefficient. 

1. Analysis of one-particle flow 

The most employed method of flow analysis at (ul- 
tra)relativistic energies is the event-plane method [231 
l24Ll34j . which relies on the event-by-event determination 
of an estimate of the reaction plane, the so-called "event 
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plane." Once the latter has been estimated (and vari- 
ous procedures to correct for acceptance issues have been 
performed) , one correlates its azimuth with that of each 
outgoing particle, assuming that all correlations between 
the event plane and a given particle (that is, actually, all 
two-particle correlations) are due to flow. Eventually one 
must correct for the event-plane dispersion (computed 
with the help of "subevents") on a statistical basis: 



v n (PT,y) = 



(cosn(-0 - vp n )) 
(cosn(* H -$ H )} : 



(Al) 



where if), ^ r, and respectively denote the azimuths of 
the particle under study, the event plane, and the real re- 
action plane, while the denominator measures the event- 
plane dispersion. 

Beside this first method, it has long been known that 
anisotropic flow can be analyzed with two-particle az- 
imuthal correlations |4jj j without having to estimate the 
reaction plane in each event. The procedure consists in 
building a two-particle correla tor, similar to that em- 
ployed in interferometry studies [10J , by forming the ratio 
of the number of "real" pairs (of particles in a same event) 
with relative angle A(f> over the number of "background" 
pairs (mixing particles from different events) with A(f>: 



N pairs {Acj)) 

-Wmixcd(A0) 



(A2) 



As usual, dividing by "mixed events" automatically cor- 
rects for acceptance anisotropies, so that one can even 
work with a detector having only limited azimuthal cov- 
erage | llj. while the event-plane method requires an al- 
most perfect azimuthal symmetry. The key point in con- 
structing C2(A0) is that its Fourier coefficients, which 
can be deduced by fitting the function, are precisely 
(cos(nA0)) = (v n ) 2 and (sin(nA^)) = 0. Letting first 
both particles in the pair run over the whole phase 
space covered by the detector, one obtains an estimate of 
"integrated flow", v n , corresponding to some (detector- 
dependent) average of the coefficient. Restricting then 
one, and only one, of the particles in the pair (whose 
azimuth will be denoted by tp) to some definite particle 
type, transverse momentum px and rapidity y, while let- 
ting the other (azimuth 0) be any particle in the event, 
one constructs a correlator whose Fourier coefficients are 



(cos n(ip - (f>)) = v n v n (p Tl y) (A3a) 
(siim(V> - </>)) = 0. (A3b) 

The second identity reflects the evenness of — <t>) 

(when parity is conserved), while the first yields the "dif- 
ferential flow" v n (pT,y)- Please note that since only 
one particle per pair belongs to a small phase space bin 
while it is correlated to all other particles in the event, 
it follows that statistical errors are the same as with the 
event-plane method. Finally, the bias from "nonflow" ef- 
fects | 42j] is of the same order of magnitude within both 



event-plane and two-particle correlation methods, but it 
is easier to subtract, when it is possible, in the latter, as 
exemplified in the case of unwanted correlations due to 
global momentum conservation in Ref. |43l ]. 

As already stated, the main limitation of both the 
event-plane and the two-particle methods is their rely- 
ing on the assumption that all azimuthal correlations 
between particles result from their correlation with the 
reaction plane 23]. In other words, they neglect non- 
flow correlations, whose magnitude is known to be large 
at ultrarelativistic energies |44j . One may try to sub- 
tract part of the nonflow effects by performing cuts in 
phase space, correlating together only particles that are 
widely separated (which certainly accounts for short- 
range correlations), but this results in larger statisti- 
cal errors, while not removing entirely all unwanted ef- 
fects. The only systematic way to remedy the problem 
of nonflow correlations in the flow analysis is to apply 
improved methods of analysis, based on multiparticle 
correlations [2^, |2^, 13(1 El E3- which have been im- 
plemented at the CERN Super Proton Synchrotron ^| 
and at RHIC |4^. The essence of these methods is that 
the relative magnitude of nonflow effects decreases, while 
that of collective anisotropic flow grows, when one consid- 
ers the cumulants of correlations between an increasing 
number of particles. Measuring cumulants of four-, six- 
particle correlations, one thus minimizes the systematic 
error due to nonflow effects, the ultimate case being the 
use of Lee- Yang zeroes |3^ , equivalent to "infinite- 
order cumulants," which isolate collective behaviors in 
the system, i.e., flow effects. The price to pay is an in- 
crease in statistical uncertainties, but the latter is mod- 
erate at RHIC and LHC energies, especially if one uses 
all detected particles in the analysis. 3 

In practical analyses, these improved methods necessi- 
tate the computation of a generating function of multi- 
particle correlations, G(z), where 2 is a complex variable 
(see Eq. (5) in Ref. 28] or Eqs. (3), (5) in Ref. 32]). One 
then derives estimates of integrated flow: in the cumu- 
lant approach, by extracting the successive derivatives 
of lnG(z) at z = and identifying them with the cor- 
responding derivatives of In Jo(2f ra |z|) [2^]; when using 
Lee- Yang zeroes, simply by looking for the location of 
the first zero of G(z) in the complex plane (see Ref. [3l| . 
Eq. (9)). Once estimates of integrated flow have been 
obtained, they are used to compute values of differential 
flow v n (pT,y) for particles in a small (pr,y) bin, whose 
azimuthal angle we shall denote by ip. This is done by 
correlating ip to the genera ting function (see Eq. (26) in 
Ref. [13 or Eq. (9) in Ref. 



3 One should not worry about possible double particle counting, 
which amount to (unphysical) nonflow effects, and are thus au- 
tomatically taken care of in the methods. 
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2. Analysis of two-particle flow 

Any of the methods recalled in the previous Section 
can also be employed to measure the anisotropic flow of 
pairs as well, modulo small modifications. Whatever the 
method, the first step is strictly the same, namely the 
construction of the event plane (and the computation of 
its statistical dispersion) in the event-plane method, or 
the measurement of estimates of integrated flow in the 
two-particle and multiparticle methods. 

We shall now describe the changes that must be made 
to the measurement of differential single-particle flow in 
order to analyze the coefficients uj? air , « p ^ r - In short, the 
first necessary modification is the obvious replacement of 
tp (the azimuth of "differential" particles) by the pair an- 
gle i?!)p a i r ; then no further change is needed to obtain the 
cosine coefficient v pair , whose measurement strictly par- 
allels that of v n (pT, y), while for 7j pa J r one should replace 
the "cos n-0" term that is correlated to either the event 
plane or the other particles or a generating function with 
a "sin n</> pair " term. Let us be more explicit: 

• In the event-plane method, the pair-flow Fourier 
coefficients are given by the averages 



( 



COS 77 



Apair 



pair 

u s,n 



(cOS77(* fl - 

(sinn(</>P air - *„)) 
(cosn(*^ - $r)) ' 



(A4a) 
(A4b) 



where is the event-plane azimuthal angle and 
the averages run over pairs and events. Note the 
analogy between Eqs. 1|A1|) and I] A4afl . 

When using two-particle correlations, one builds a 
two-point correlator C2(0 pair — <j>), where <j> is any 
particle in the same event as those involved in the 
pair, with the trivial exception of the pair parti- 
cles to avoid autocorrelations. In opposition to 
the correlator used in single-particle flow studies, 
Ci (</> pair — cf>) is no longer an even function, so that 
its Fourier expansion has nonvanishing both cosine 
and sine coefficients which are related to the pair- 
flow coefficients, namely 



(cOS77((/) Pair -(f))) = 



(si 



sin 77 



/pair 



)) 



u n v s,n J 



(A5a) 
(A5b) 



where v n is the integrated flow while the aver- 
ages run over all (0 pair , <f>) in each event, then over 
events. Once again, Eq. (|A5a|) is reminiscent of 
Eq. i|A3a|) . while the difference between Eqs. (|A3bp 
and i|A5b|l is due to the fact that whereas single- 
particle emission is symmetric with respect to the 
reaction plane, pair emission, on the other hand, is 
not symmetric, see Sec. Inland Fig. ^ 

To measure the Fourier coefficient w pa ^ r in the Lee- 
Yang zeroes method, one should replace cos n(ip—9) 



by sinn(0 P air — 9) m the numerator of Eq. (12) 
(resp. Eq. (9)) in Ref. P (resp. Ref. P) 4 

• Finally, in the cumulant method, the relevant cu- 
mulants when measuring u p£ JJ r are the imaginary 
parts in the power-series expansion of Eq. (26-27) 
in Ref. |2jj , while the real parts are needed for the 
analysis of single-particle differential flow v n (px, y) 
or of the cosine coefficients v pa ". As a result, the 
interpolation formula that allows one to extract 
the cumulants is similar to Eq. (B7) of Ref. [2^ 
(resp. Eq. (11) of Ref. 29]), modulo the replace- 
ment (X Piq ,Y Piq ) — > (Y Piq ,—X Piq ), where X p>q and 
Y p>q are still given by Eq. (B6) in Ref. HJ. 

Even if every method of single-particle flow analy- 
sis can in principle also be used to measure pair flow, 
modulo the modifications we described above, there ex- 
ists a clear difference between the two-particle methods 
(both event-plane and two-particle correlation methods) 
on the one hand, and the multiparticle approaches on 
the other hand. As a matter of fact, we already men- 
tioned that both two-particle methods rely on the as- 
sumption that all correlations between two arbitrary par- 
ticles are due to anisotropic flow; in other words, that 
other sources of two-particle correlations are absent, or 
at most weak j^. Now, if the purpose of using one of 
these methods is precisely to measure some azimuthally- 
dependent two-particle effect, the procedure is somehow 
self-contradictory. Such an inconsistency does not af- 
fect the measurement of pair flow through multiparticle 
methods, since the latter do not assume that two-particle 
correlations are inexistent, they merely minimize their ef- 
fect. 



APPENDIX B: SYMMETRY PROPERTIES OF 
PAIR-FLOW COEFFICIENTS 

In this Appendix, we list a few mathematical prop- 
erties of the pair-flow coefficients 7j pair for the sake of 
completeness. 

The invariance of the two-particle distribution under 
the transformation (</>i,</>2) — ► (4>i + 2tt, 02 ) translates 
into the (pseudo) periodicity property 



pair 



(A0 + 2vr) 



pair 



(A0)e 



-2i7rnx 



(Bl) 



where x has been defined in Eq. If x is changed 



to x' in Eq. (3J|, 



pair 



(A0)e 



pair 



\n{x — x )A0 



(A</>) is changed to v' n pa " (A0) 



If the system has symmetry with respect to the re- 
action plane (no parity violation), the two-particle dis- 
tribution is unchanged under the joint transformation 



4 We assume for simplicity that m takes the value 1 in the cited 
equations. Values m > 1 correspond to higher harmonics, for 
which the same modifications apply. 
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, Acj>) -> {-<tP ai \ -A0). At the level of the v p s f{-A(f>) = v^ n (Acj)) = -t^(A^). (B3) 

Fourier coefficients, this symmetry gives v? 3 -" (~ A<fi) = 
v^'(A</)). Together with property this yields 

< air (-A0) = ^7(A4>) = K ah '(A0)]*. (B2) 

The corresponding properties for the real Fourier coeffi- 
cients v^" and v^ a " are These various properties may prove useful to check that 

measured estimates of the Fourier coefficients behave 
^ r (-A0) = v^ a i T n (A(/)) = v^ T (A(f>) "properly". 
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